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Abstract 

We propose a novel approach for efficient tuning of optical properties of a high refrac¬ 
tive index subwavelength nanoparticie with a magnetic Mie-type resonance by means of fem¬ 
tosecond laser irradiation. This concept is based on ultrafast photo-injection of dense (> 10^° 
cm^^) electron-hole plasma within such nanoparticie, drastically changing its transient di¬ 
electric permittivity. This allows to manipulate by both electric and magnetic nanoparticie 
responses, resulting in dramatic changes of its scattering diagram and scattering cross section. 

We experimentally demonstrate 20% tuning of reflectance of a single silicon nanoparticie by 
femtosecond laser pulses with wavelength in the vicinity of the magnetic dipole resonance. 
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Such single-particle nanodevice enables to design fast and ultracompact optical switchers and 
modulators. 

All-dielectric "magnetic light" nanophotonics based on nanoparticles of high refractive index 
materials allows manipulation of a magnetic component of light at nanoscale without high dissi¬ 
pative losses, inherent for metallic nanostructures.^^ This "magnetic light" concept has been im¬ 
plemented for nanoantennasphotonic topological insulators,!^ broadband perfect reflectors,!^ 
waveguides,!^ cloacking,!^ harmonics generation,!^ wave-front engineering and dispersion con¬ 
trol.!^ 

Such magnetic optical response originates from circular displacement currents excited inside 
the nanoparticle by incident light. This opens the possibility of interference between magnetic and 
electric modes inside the dielectric nanoparticle at certain wavelength. One of the most remark¬ 
able effects based on this concept is formation of the so-called Huygens source, scattering forward 
the whole energy ,!i^ while for another wavelength range the nanoparticle can scatter incident light 
almost completely in backward direction.! ^^!^^ ! Therefore, manipulation by both electric and mag¬ 
netic resonances paves the way for efficient tuning of the dielectric nanoparticle scattering in the 
optical range. The spectral positions of the electric and magnetic dipole resonances depend on 
the particle shape and environment.! ^ ! ^ ! ^ ! ^^!^^!^^ ! Alternatively, these resonances can be tuned per¬ 
manently via changes in dielectric permittivity, which was achieved by annealing of amorphous 
silicon nanoparticles .!^ 

However, modern optical technologies require fast, large, and reversible modulation of optical 
response of ultracompact functional nanostructures. For this purpose, different types of optical 
nonlinearities both in metallic!^ and dielectric structures!^ were utilized such as Kerr-type nonfln- 
earities,!^^!!^free carriers generation!^2Hl2!and variation of their temperature,!^ as well as relatively 
slow thermal nonlinearity.!^ Since plasmonic structures have high inherent losses, while photonic 
crystals or graphene-based structures!^ are much larger than the wavelength, it is advantageous 
to use low-loss and subwavelength high-index particles with electric and magnetic responses (the 
"magnetic light" concept). Moreover, nonlinear manipulation by scattering characteristics (power 
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Figure 1 : A schematic illustration of scattering manipulation by a single weak (a) and intense (b) 
femtosecond laser pulse. Intense femtosecond laser pulse switch-on a Huygens source regime, 
when incident light is scattered in the forward direction. 

pattern and cross section) of a nanostructure via its magnetic response gives an additional efficient 
tool for ultrafast all-optical switching and routing. 

In this work we propose a novel approach for manipulation by the electric and magnetic re¬ 
sponses of a high-index dielectric nanoparticle, employing its ultrafast photoexcitation by fem¬ 
tosecond laser irradiation. Specifically, we demonstrate theoretically possibility of large tuning of 
scattering properties of the single nanoparticle (near the regime of the Huygens source as shown 
in Fig. and achieve experimentally 20 % changes of reflection from silica surface with a silicon 
nanoparticle under femtosecond laser irradiation. 

Femtosecond (fs) laser pulses are known to provide strong photo-induced electronic excitation 
of free carriers in diverse materials (owing to usually low electronic heat capacities), which is ac¬ 
companied by sub-10 fs thermalization of free carriers and dramatic variation of material optical 
characteristics during the pumping fs-laser pulse (usually < 100 fs). Such almost prompt, fs-laser 
induced optical tunability appears to be much broader for (semi)insulating materials with very 
minor initial carrier concentrations, extending in a sub-ablative regime in the visible and near-IR 
ranges. In detail, such ultrafast transient modulation of optical dielectric permittivity in semicon¬ 
ductors and dielectrics is related to transient variation of free-carrier (electron-hole plasma, EHP) 
density peh through its basic intraband and interband contributions Simultaneously, ultrafast 
transient optical modulation is additionally enhanced due to a strong prompt EHP-driven isotropic 
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renormalization of direct bandgap, resulting in drastic enhancement of interband transitions and 
corresponding red spectral shift of the optical dielectric permittivity Remarkably, such strong 
optical modulation in semiconductors, requiring high EHP density peh > 10^® cm is rapidly 
reversible on a picosecond time scale, e.g., down to ~6-7 ps for c-Si owing to three-body Auger 
recombinationl^^EQl or even down to ~1 ps for a-Si ® 

In this work a realistic pump-pulse averaged dependence of the dielectric permittivity for the 
photo-excited silicon versus incident fs-laser fluence at 800-nm laser wavelength was obtained 
through modeling and fitting of experimental data from our previous worlJ^ (Fig. [^, which is 
well consistent with data from many other papersparticular, we analyze single-shot fs- 
laser pump self-reflectivity (R) from an atomically smooth silicon surface at its s- (7?j(45°)) and p- 
{Rp{45°)) polarizations at the 45°-incidence angle and variable effective (absorbed) laser fluences 
^eff = (1 ~^s,p(45°,F))-F, where F is the incident fluence, using a model transient dielectric per¬ 
mittivity for photo-excited silicon. Commonly, such model dielectric permittivity, being a function 
of peh, can be expressed as a sum of interband- and intraband-transition based terms 


e(®,Peh) 


£ib(C0*) 


Peh \_ ft>pi(Peh) /_ i 

pbf/ (0^+l/{x^{psh))\ (OXeipeh) 


( 1 ) 


where the above mentioned peh-dependent bandgap shrinkage effect on interband transitions is 
accounted by introducing effective photon frequency CO* = to-f 0peh/Pbgr with the factor 0, the 
characteristic renormalization EHP density Pbgr ~ 1 x 10^^ cm^^,^ being typically about 5% of 
the total valence electron density 2x10^^ cm~^ in Si) to provide the ultimate 50% electronic 
direct bandgap renormalization,!^ i.e., Ti& ^ IJ eV of the effective minimal gap 3.4 eV in 
silicon,!42E4l^]^qe is characteristic band capacity of the specific photo-excited regions of the 
first Brillouine zone in the k-space (e.g., pbfiL) ~ 4 x 10^^ cm^^ for E-valleys and pbf(2f) ~ 4.5 x 
10^^ cm^^ for X-valleys in Si), affecting interband transitions via the band-filling effect .1^2^^221391 
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Figure 2: Derived changes of real (black curve) and imaginary (red curve) components of optical 
dielectric permittivity versus Fgff at the 800-nm pump wavelength relatively their initial values. 
Inset: experimental pump self-reflectivity dependencies i?s,p(45°,Feff) (red and blue dots, respec¬ 
tively, adapted fromESi) on effective laser fluence with their corresponding fitting red and blue 
model curves i?s,p(45°,Peh)- The black dashed line shows a typical value of laser fluence, which is 
enough for strong changing of optical properties of bulk silicon or a silicon nanoparticle (see Fig. 
4) at sub-damage conditions. 


The bulk EHP frequency cOpi is defined as 


£oe|,f(Pel.)m*p,(Peh) ’ 

where the effective optical (e-h pair) mass m*pj 0.14me in L-valleys or 0.19 in X-valleys^^MlIllIll] 
is a peh-dependent quantity, varying versus transient band filling due to the band dispersion and 
versus bandgap renormalization.^SThe high-frequency electronic dielectric constant e^f was mod¬ 
eled in the form £hf(peh) = 1 +£hf(0) x exp(—Peh/Pscr), where the screening density Pscr ~ 1 x 
10^^ cm^^ was chosen to provide £hf —)■ 1 in dense EHP. The electronic damping time tg in the 
regime of dense EHP at the probe frequency (Opr was taken, similarly to metals, in the random 
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phase approximation as proportional to the inverse bulk EHP frequency COpj 


( 1284 ^ 

1 -Fexp 

hco 

ksTe 

\;r2v^(0pi ) 

{nksTeY + {fi(o) 


where Ep ~ 1 - 2 eV is the effective Fermi-level quasi-energy for electrons and holes at peh < 
1 X 10^^ cm^^, h and kg are the reduced Planck and Boltzmann constants, respectively, and 
is the unified EHP temperature, being a weak function of Peh® Here, the latter relationship was 
evaluated for tico > in the form Te(peh) ~ 3 x 10^/((Opi(peh)), accounting multiple carrier scat¬ 
tering paths for the three top valence sub-bands, and multiple X-valleys in the lowest conduction 
band of silicon. 

The resulting Peh-dependent oblique-incidence pump reflectivities i?s,p(45°,peh) model, calcu¬ 
lated using common Fresnel formulae, fit well the extracted experimental reflectivity dependences 
i?s,p(45°,Feff) in Fig. 2 with the characteristic initial dip and the following rise. Such reasonable 
fitting in Fig. [^provides an important relationship between magnitudes peh and Feff in the region, 
covering the reflectivity dip and rise, which were used to plot the derived optical dielectric permit¬ 
tivity versus Feff (Fig. |^. 

Such large changes of optical dielectric permittivity properties in silicon at fluences below its 
melting (Fgff ~ 0.17 J/cm^l^ and ablation thresholds (for spallation under these experimental 
conditions, Feff ~ 0.3 J/cm^l^ can significantly alter optical response of a silicon nanoparticle, 
supporting a magnetic Mie-type resonance. Using the extracted dielectric permittivity values of 
photoexcited silicon, we study such optical tuning of scattering properties of a silicon nanoparticle 
in vacuum near its magnetic resonance by means of full-wave numerical simulations carried out in 
CST Microwave Studio. 

We numerically analyzed optical properties of a spherical (the diameter D = 2r = 210 nm) 
silicon particle with its dilectric permittivity, depending on laser fluence as shown in Fig.[^. The 
chosen nanoparticle diameter corresponds to excitation of a magnetic dipole Mie-type resonance 
in the vicinity of the femtosecond laser wavelength A ~ 800 nm. Its scattering cross-section and 
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Figure 3: Numerically calculated scattering spectrum of a silicon sphere of a diameter D = 210 nm 
(a). Insets: numerically calculated electric and magnetic fields distributions inside and near the 
silicon sphere. Scattering diagrams of the silicon sphere at A = 740 nm (b), 800 nm (c) and 
860 nm (d). The black arrows in (b-d) indicate the direction of light incidence. 


scattering diagram are well-known to be rather spectrally sensitive in the vicinity of the magnetic 
resonance.12^ In particular, at some wavelengths, where magnetic and electric dipoles induced in 
the nanoparticle are almost equal and oscillate in phase, the silicon nanoparticle works as a Huy¬ 
gens source with suppressed backward scattering.1211^ On the other hand, the scattering diagram 
can be tuned to the regime of suppressed forward scattering ("reverse" Huygens source), when 
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Figure 4: Numerically calculated normalized scattering spectra (a) and scattering diagram at A = 
800 nm (b) of the 210-nm silicon particle irradiated ati^ff = 0 (black curve), 35 mJ/cm^ (red curve), 
and 100 mJ/cm^ (blue curve). The dashed line corresponds to the wavelength of the pumping 
femtosecond laser. 

magnetic and electric dipoles oscillate with the phase difference of TtH^^ 

The results of numerical simulations of scattering and scattering properties of the silicon nanopar¬ 
ticle in the vicinity of the magnetic dipole resonance are shown in Fig. The calculated distribu¬ 
tions of electric and magnetic fields indicate an emerging magnetic dipole moment in this range 
(Fig.|^), while electric dipole moment is weak in the vicinity of 800 nm. At the points "b" and "d" 
the internal field represents the fields superposition of electric and magnetic nonresonant dipole 
modes. Since both scattering spectra (Fig. [^) and scattering diagram (Fig. |^-d) are strongly 
wavelength-dependent, it is possible to tune these two parameters by varying its dielectric permit¬ 
tivity. 

Significant tuning of optical properties in the silicon nanoparticle near its magnetic dipole reso¬ 
nance requires time for the internal electromagnetic mode formation shorter than the corresponding 
electronic damping (thermalization) time and laser pulse duration. For the 210-nm silicon nanopar¬ 
ticle, the magnetic resonance mode has approximately the 50-nm full-width at the half- maximum, 
and, therefore, the Q-factor is about 16, corresponding to the mode formation time in such open 
resonator of about 7 fs. Specifically, for dense EHP with (peh > 10^^ cm ^ thermalization of free 
carriers proceeds over the characteristic plasma density-dependent times Tg ~ 10-100 fs (see Eq. 
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(3)) during the 100-fs pumping laser pulse. Henee, ultrafast plasma density-dependent ehanges in 
optieal dieleetrie permittivity of the partiele eontinuously oeeur during its laser pumping and are 
almost adiabatieally followed by the internal eleetromagnetie mode. 

To simulate numerieally ehanges in optieal properties of the photoexeited silieon sphere with 
the diameter D = 210 nm at the wavelength A = 800 nm, we use the derived above dependeneies 
of ARe(e) and Alm(e) on absorbed laser fluenee Fgff (Fig.[^). The eonsidered range of absorbed 
fluenees F^ff <100 mJ/em^ eorresponds to non-destruetive regime of the laser-partiele interaetion, 
but laser fluenees are still high enough to generate rather dense EHP (peh~lx 10^^ em for 
effleient switehing of the nanopartiele optieal properties. As was mentioned above, the seattering 
diagram of the 210-nm nanopartiele is almost symmetrie at A = 800 nm and Feff ~ 0 (Fig. |^), 
while the seattering eross seetion (Cext) normalized on has rather high value of about 9 (see 
Figs. and|^). The latter parameter is ehanged almost by three times with fluenee inereasing 
up to Feff =100 mJ/em^ at the fixed wavelength of 800 nm, owing to the strong shift of the peak 
position of the seattering speetrum (Fig.|^). Its seattering diagram appears to be also very sensitive 
to the eorresponding ehanges of the dieleetrie permittivity. At A ~ 800 nm, the transition from 
the typieal dipole seattering diagram to the Huygens souree forward seattering is observed in the 
fluenee range Feff = 0—100 mJ/em^ (Fig. |^). 

In order to prove the eoneept of manipulation by the eleetrie and magnetie responses of a 
nanopartiele via optieal photo-injeetion, we performed experiments with silieon nanopartieles of 
different sizes. First, we have fabrieated arrays of truneated eonieal silieon nanopartieles with 
period A 800 nm from an silieon film with thiekness h=220 nm deposited by PECVD on a fused 
silica substrate by means of inductive coupled plasma etching through electron-beam lithography 
prepared metal mask. The base (r^) and top (r^) radii of the nanopartieles were designed basing on 
preliminary numerical modeling, predicting the spectral position of the magnetic dipole resonance. 
Two different arrays were studied, composed by "near-resonance" (r^, 120 nm and 70 nm) 

and "off-resonance" {rj, 70 nm and Vt ~ 20 nm) nanopartieles. It should be noted here that the 
condition of Huygens source has been observed previously in visible and IR ranges for spherical. 
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cylindrical and even eonieal semiconduetor (silieon, GaAs, etc.) nanoparticles 

The broadband speetral measurements of the refleeted signal from a single nanopartiele (Rs) 
were earned out by means of strong focusing and collection of light (A =400-900 nm) to a spec¬ 
trometer (Horiba LabRam HR) through an aehromatic objeetive with a numerieal aperture NA=0.95. 
This objeetive allows to irradiate and eolleet light from an area of diameter ~ 1.22X/NA < A. 
The resulting refleetion spectra from different nanoparticles exhibit pronounced maxima (reso- 
nanees) in different speetral ranges: one resonance is near A = 500 nm for the near-resonanee 
nanoparticle and two resonanees at A = 700 nm and A = 850 nm for the off-resonance nanopartiele 
(Fig.|^). Our numerieal modeling of refleetion from silieon nanopartieles on Si02 substrate with 
the given dimensions and tabulated speetral dispersions^^ reveals, that the observed maxima for 
near-resonance nanopartiele eorrespond to eleetrie (700 nm) and magnetic (810 nm) dipole res¬ 
onances, while the maximum in off-resonanee nanopartiele spectrum has the magnetie response 
origin (Fig.[^). 

The nonlinear measurements of fluenee dependent ehanges in the refleetion signal from a single 
nanopartiele were earried out by using the objeetive (NA=0.95), similarly to the linear measure¬ 
ments. A eommereial femtoseeond laser system (TiF-lOOF, Avesta Pojeet) was used as an intense 
light source, providing 100-fs laser pulses at 790 nm central wavelength (FWHM 20 nm), with 
maximum pulse energy of 5 nJ and pulse duration of 100 fs at the repetition rate of 80 MHz. 
Laser energy was varied and eontrolled by an aeousto-optieal modulator (R23080-3-LTD, Gooch 
and Housego) and a power meter (FielfMax II, Coherent), respeetively. The ehosen energy range 
eorresponds to fluenees less than 100 mJ/em^, whieh is elose but well below the melting threshold 
(F^ 250 mJ/em^l^ and the observed damage fluenee (F ^ 100 mJ/em^) of the nanoparticles. 
The refleeted femtosecond laser pulses were colleeted by the focusing objective and directed to a 
ealibrated CCD eamera to measure refleeted average fluenee with high spatial resolution. 

The wavelength of the femtoseeond laser irradiation is loeated at the blue shoulder of the mag¬ 
netie dipole resonanee of the near-resonanee nanopartiele, while it is far away from all resonanees 
of off-resonance nanopartiele (Fig. |^,b). The comparison of relative variation of reflectivity from 


10 



Figure 5: Experimental (a) and numerical (b) reflection spectra of a truncated conical silicon 
nanoparticle (see the bottom left inset) with radii 70 nm and 20 nm (green curve), or 
~ 120 nm and 70 nm (red curve) with h = 220 nm; spectrum of the femtosecond laser 
(right axis). The upper left insets in (a): optical image of silicon nanoparticles arrays. The bottom 
left inset in (a): SEM image of the "near-resonance" nanoparticle (the scale bar is 50 nm). The 
red filled curve in (a) represents the experimental spectrum of the femtosecond laser. The insets 
in (b) represent calculated electric field distributions near and inside the nanoparticles with the 
abbreviations: MD - magnetic dipole resonance, ED - electric dipole resonance, (c) Experimental 
(dots) and theoretical (solid lines) dependencies of normalized reflectance change on laser fulence 
for 220-nm thick silicon film (black), the "near-resonance nanoparticle" (red) and "off-resonance 
nanoparticle" (green). The values of Rso are different for each sample and correspond to E = 1 
mJ/cm^. The upper left inset: a schematic illustration of reflectance measurements from a single 
nanoparticle 
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the surface with the near-resonance nanoparticle, the non-resonant nanoparticle, and bulk silicon 
on laser fluence reveals the strongest response from the near-resonance nanoparticle. Specifically, 
the reflectance of the near-resonance nanoparticle exhibits 20% relative growth with incident flu¬ 
ence increasing up to ~ 30 mJ/cm^, while the reflection signal from the 220-nm thick silicon 
film goes down and the off-resonance nanoparticle reflection demonstrates intermediate behavior 
(Fig.ig^). 

Such changes in reflection dependencies on laser fluence are totally governed by photo-injection 
of EHP, because we found good agreement of experimental dependencies of ARs/Rso on incident 
fluence with the numerical simulations of changes in reflection signal (Fig.[^), taking into account 
the calculated dielectric function dependence on laser fluence (Fig.[^. The growth of the reflectiv¬ 
ity of the near-resonance nanoparticle is caused by the blue shift of the magnetic dipole resonance 
owing to the generation of FHP, decreasing dielectric permittivity of silicon with increase of laser 
fluence (Fig.[^. 

In summary, photo-injection of dense electron-hole plasma in a dielectric nanoparticle, sup¬ 
porting a magnetic dipole resonance in the optical range, paves the way for effective light manip¬ 
ulation on subwavelength scale not only by tuning of scattering cross section of the nanoparticle, 
but also by tuning of its scattering diagram. In the frame of this concept, the 20% switching of 
reflection from a silica surface with a silicon nanoparticle photoexcited by a femtosecond laser 
pulse has been shown, enabling high-efficient light manipulation on the subwavelength scale. In 
general, since the transient electronic dynamics under VIS-IR femtosecond laser photoexcitation is 
governed by similar physical processes (FHP generation, bands filling, bandgap renormalization, 
and ion potential screening) in different semiconductors, the electron-hole plasma induced tuning 
of magnetic resonance can be applicable for broad range of materials and wavelengths. Therefore, 
tuning of dielectric permittivity of a single dielectric nanoparticle via photoexcitation of dense 
electron-hole plasma opens a novel class of ultracompact nanodevices with potentially ultrafast 
timescale, based on the diversity of "magnetic light" effects. 
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